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Abstract 

We describe an efficient approximation algorithm for evaluating the ground-state en- 
ergy of the classical Ising Hamiltonian with linear terms on an arbitrary planar graph. 
The running time of the algorithm grows linearly with the number of spins and expo- 
nentially with 1/e, where e is the worst-case relative error. This result contrasts the well 
known fact that exact computation of the ground-state energy for the two-dimensional 
Ising spin glass model is NP-hard. We also present a classical approximation algorithm 
for the Local Hamiltonian Problem or Quantum Ising Spin Glass problem on a planar 
graph with bounded degree which is known to be a QMA-complete problem. Using a 
different technique we find a classical approximation algorithm for the quantum Ising 
spin glass problem on the simplest planar graph with unbounded degree, the star graph. 

1 Introduction 

Ising spin glasses model physical spin systems with random, competing interactions due to 
disorder in the system [II. In order to make meaningful predictions about such systems 
one can consider statistical ensembles that represent different realizations of the couplings. 
For a particular realization of the couplings one is generally interested in finding algorithms 
to determine properties such as the spectrum, partition function, the ground-state or the 
ground-state energy. Algorithms with a running time that is a polynomial in the problem 
size are called efficient, in contrast with inefficient procedures that take super-polynomial or 
exponential running times. Connections between disordered systems in statistical physics and 
questions of computational complexity have been extensively explored, see e.g. [2] and [3]. 

In this paper we find a new application of computational complexity tools to spin glass 
problems, namely a rigorous approximation algorithm to determine the ground-state energy of 
a classical or quantum Ising spin glass on a planar graph. It has been shown that to determine 
the ground-state energy of a classical Ising spin glass on a 2D lattice with linear terms exactly 
is computationally hard, or NP-complete [4]. Terhal and Oliveira showed that determining 
the smallest eigenvalue of a quantum Ising spin glass of n qubits on a planar graph with 
l/poly(n) accuracy is QMA-complete [5]. Aharonov et al. [S] showed that even determining 
the smallest eigenvalue for audits on a line is QMA-complete. Therefore the running time of 
any algorithm for computing the ground state energy exactly must scale super-polynomially 
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with the number of spins n (under the assumption Py^NP). Thus it is natural to look for 
approximation algorithms that solve the problem in polynomial time at the cost of providing 
a slightly non-optimal solution. 

Let us define the (quantum) Ising spin glass problem precisely. Let G — (V, E) denote any 
graph with n vertices and let u, v € V be vertices of the graph. One can associate a 2-local 
Ising spin glass Hamiltonian H with interaction graph G — (V, E), 

h= 53 Qu,v + J2 L -» = Q+L - w 

(u,v)eE u£V 

Here Q u ,v is quadratic in Pauli operators and L u is linear. We assume that ||<2 U) „||, \ \L U \\ < 
poly(n). It is also assumed that Ty H — 0. We obtain the classical Ising spin glass, the 
function H(S), by setting Q UjV = c uv S u S v and L u = d u S u for classical spins S u = { — 1, 1}, 
i.e. 

H(S) = 2 c uv S u S v + ^ d u S u (2) 

u,v u 

In the classical case c uv and d u can be given as some m-bit numbers. Bieche et al. [7] has 
shown that for planar graphs G the problem of determining the minimum value of H(S) (the 
ground-state energy) and the associated assignment S can be solved efficiently if there are no 
linear terms (i.e. all d u = 0). 

Our approximation algorithm for the classical and quantum Ising spin glass is relevant 
in light of the recent research on quantum adiabatic approaches for finding the ground-state 
of a classical or quantum Ising spin glass and solving other NP-complete problems. The 
paradigm of adiabatic quantum computation (AQC) was first introduced in [8]. An adiabatic 
computation proceeds by slowly varying the system's Hamiltonian starting from some simple 
Hamiltonian Hq at the time t = and arriving to, for example, a classical Ising spin glass 
Hamiltonian H at t = T (the final Hamiltonian can also capture other NP-hard problems). 
The adiabatic theorem, see e.g. [9] , guarantees that if one starts from the ground state of Hq 
and the running time T is large compared to the inverse spectral gap at all times then the 
final state is close the ground-state of H. 

If we assume the validity of the conjecture that a quantum computer can not solve NP- 
complete nor QMA-complete problems (see for classical spin glasses the analysis in [lOllllj ) it 
follows that AQC is a means to obtain an approximation to the ground-state and the ground- 
state energy. The quality of this approximation and its dependence on the hardness of the 
problem and the adiabatic path are at present not well understood. Some experiments on 
a physical realization of the classical Ising model in a transverse field show that quantum 
annealing can lead to faster equilibration of the system [12]. Also, it was recently shown that 
a quantum algorithm can provide a square-root speed-up over a classical simulated annealing 
algorithm [15] . These results do not show that the Ising spin glass problem becomes easy 
on a quantum computer, but point to possible advantages of using a quantum computer to 
obtain efficient approximation algorithms. Recently the company D-wave has claimed to have 
implemented the Ising spin glass Hamiltonian (with additional edges on the diagonals of the 
lattice) and an adiabatic evolution for 16 qubits on a 4 x 4 square lattice, see [T3]. The hope of 
this endeavor is that such system outperforms classical computers in (approximately) solving 
optimization problems. 

Given these claims and results about quantum speed-ups it is clearly interesting to consider 
how well an approximation to the ground-state energy can be obtained by purely classical 
algorithmic means. 
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The area of approximation algorithms is an active area of research in computer science, 
see e.g. [15] . Such algorithms are often of practical importance for generic hard problems 
for which we are willing to trade off the relative quality of the approximation versus the 
running time of the algorithm. Three main types of approximations to optimization problems 
can be distinguished. A problem is said to have a polynomial time approximation scheme 
(PTAS) if given any e > 0, there is an algorithm A e which for any instance I produces a 
solution within (1 ± e) times the optimal solution. In addition A e has a running time which 
is a polynomial in the input size of /. Observe that the running time of a PTAS is only 
required to be polynomial in input size, and it can have an arbitrary dependence on e (for 
example n 1 ^ or poly(n) are valid running times for a PTAS). A stronger notion is that 
of a fully polynomial time approximation scheme (FPTAS), where the running time of the 
approximation scheme is required to be polynomial both in the size of the input and in (1/e). 
For the classical Ising spin glass problem an instance is a particular graph and set of values 
of the weights c uv and d u . The optimal solution is the minimum value of the energy H(S). 

For the Ising spin glass problem on planar graphs, we can exclude the possibility of a 
FPTAS (assuming P ^ NP). This is due to the fact that the problem is NP-hard even when 
c uv , d u are restricted to be either —1,0 or +1 4 ■ Let us assume that we have a FPTAS and 
set e = d/poly(n) for some constant 5. This gives a polynomial-time algorithm to approximate 
miris H(S) with an error which is at most 5mins H(S)/j>o\y(n). This is sufficient accuracy 
to solve the NP-complete problem exactly since H(S) is at most poly(n) and the difference 
between the minimum of H(S) and the value right above it (i.e. the energy gap) is at least 
1, that is, independent of n. 

A class of classical optimization problems for which the objective function can be efficiently 
approximated with a relative error e for some fixed e is called APX. It is known that some 
problems in APX do not have a PTAS (under the assumption P^NP), for example, Minimum 
Vertex Cover and Maximum Cut problems, see [16 . In other words, for such problems no 
polynomial-time algorithm can make the relative error smaller than some constant threshold 
value eo- Problems having this property are called APX- hard. One can use the relation 
between the Ising spin glass problem and the Maximum Cut problem, see [17], to show that 
the former is APX-hard if defined on arbitrary graphs. This is the reason why the present 
paper focuses only on planar graphs. 

Note that the existence of a PTAS for a Hamiltonian H does not imply the existence 
of a PTAS for the trivially related problem of finding the ground-state energy of H + al, 
where a is a constant; this is because the PTAS produces a solution with small relative error. 
Our PTAS for the quantum and classical Ising spin glass problem can be translated to an 
approximation algorithm with an absolute error at most eW, where W = J2^ u v \ ||Quu|| ( see 
e.g. Eq. ([6]) for the classical error analysis). Clearly, while comparing the quality of an 
approximation obtained using AQC and the classical PTAS the relevant figure of merit must 
be an absolute error, because the Hamiltonians H and H + al are physically equivalent. Note 
that when we state our running times for the approximation algorithms, we state the worst- 
case running time depending on some guaranteed error-bound. In practice, running times 
may be much faster if heuristic methods are used within the approximation algorithm. Such 
heuristic methods cannot guarantee an error-bound, but may work well for average-case or 
'real-life' instances. 

We will first consider the classical Ising spin glass on a graph G which is a two-dimensional 
lattice, and give a PTAS for this case. It has a running time T = 0(n4~), see Section [2Tl The 
intuitive idea behind this construction is simple. Assume, for simplicity, that all couplings 
between spins have similar strength. Then one can subdivide a 2D lattice into blocks of size 
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L x L by omitting the edges connecting these subblocks. The total contribution of these 
omitted boundary edges scales as 4L x — 0(n/L) and hence for large, but constant, 
L = 1/e the error that one makes by omitting these edges is bounded by at most O(en). This 
proves that there exists an approximation algorithm with absolute error. However one can 
show that the ground-state energy scales with n (see e.g. the rigorous Theorem [J) and thus 
the error will be proportional to the ground-state energy which is exactly what is desired for 
the PTAS. 

In our formulation of the problem, not all edges on the 2D lattice have similar strength, 
hence somewhat more elaborate arguments are needed to show the existence of a PTAS. 

For general planar graph this situation is more involved. Vertices in the graph can have 
arbitrary high degree and it is not clear how to divide up the graph into sub-blocks with small 
boundaries. Let W — J2( u ,v) \ \Qu,v\\ f° r the quantum Isin spin glass and W = Yl( u ,v) l c ™l 
for the classical Isin spin glass, see the definitions of the Hamiltonians in Eqs. (|1I2[) . In the 
'classical' Theorem[T]and its quantum counterpart, Theorem[2J we will show that the ground- 
state energy of an Ising spin glass on a planar graph is less than —cW for some constant 
c. This rigorously expresses the intuitive physical notion that the ground-state energy is 
extensive in the system size n. The idea of the PTAS is then as follows. We take out a 
subset of edges in the planar graph for which (i) one can show that they contribute at most 
eW to the Hamiltonian and (ii) by taking out these edges one ends up with a set of simpler 
disconnected graphs on which one can solve the problem efficiently (in poly(n)2°^ 1 / e ) time). 
Let H denote the (classical or quantum) Hamiltonian for the problem where we have taken out 
these edges. Since we only take out edges (and no vertices), a state with minimum energy for 
H is also a state for H . Let X(H) denote the ground-state energy of H. By Weyl's inequality 
\X(H) - X(H)\ < \\H - H\\ < eW. Using Theorems [1] and [2] we can relate W back to the 
lowest eigenvalue of H and hence show that the PTAS outputs a (quantum) state which has 
energy at most 0(eX(H)) higher than the true ground-state energy. 

How do we take out edges from the original planar graph? In the classical case, see 
Section l2~5l we take out edges related to a so-called outerplanar decomposition of a graph. In 
this way we end up with disconnected graphs which have a constant tree-width. It is known 
how to solve the classical spin problem on graphs with bounded tree-width (using dynamic 
programming). 

In the quantum case we cannot chose this procedure since the quantum problem on a 
graph with bounded-tree width, or even on a tree, can still be hard, see [5J. This points to 
an interesting difference between the quantum and the classical Isin spin glass. 

If, in the quantum case, we restrict ourselves to graphs with bounded-degree, we can 
apply a procedure that removes edges and leaves a set of disconnected graphs each of which 
has constant size (related to e) , see Section 13.11 Determining the ground-state energy of a 
Hamiltonian in a space of constant dimension can be done classically. Note that this PTAS 
outputs a classical description of a quantum state which has an energy 0(eA(ff))-close to the 
true ground-state energy. 

Our last result finds a PTAS for the quantum problem on a star graph, see Section [5] 
The complexity of the local Hamiltonian problem on this graph is not known. For the star 
graph it is not clear how to take out edges without introducing a large error. Hence we will 
use a different technique which uses symmetry and the rounding of interactions. In effect, we 
construct a Hamiltonian H for which X(H) and its ground-state can be determined efficiently 
and \\H — H\\ < eW. The construction works as long as all terms in H have norms in a 
range [a, 1] for a constant a. Such condition was not present in the other PTAS constructions. 
Extensions of this technique may be important for addressing the general quantum problem 
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on planar graphs. 

We note that our technique for the classical planar graph problem, i.e. using an outerplanar 
decomposition of the graph, is fairly standard for solving hard classical problems on planar 
graphs. In fact, many problems admit a PTAS on planar graphs even though approximating 
them on general graphs is known to be NP-hard, see [TH [19l [20] . While our techniques are 
similar to those of [19] and [20] at a conceptual level, our results do not follow directly from 
their work and require some new ideas. The main difficulty is that the Hamiltonian involves 
both positive and negative terms which can possibly cancel out. 

2 Classical Hamiltonians on planar graphs 
2.1 The 2D Lattice Case 

Consider the classical Ising spin glass Hamiltonian Eq. @ defined on a 2D square lattice of 
size yfn x y/n. Let e > be a fixed small constant. Without loss of generality let us assume 
that t = 1/e is an integer. For i = 0, . . . ,t — 1, let Xi denote the set of vertices u = (x, y) on 
the horizontal lines defined by {y = i (modulo t) }. Similarly, let Yi denote the set of vertices 
on the vertical lines defined by { x = i (modulo t) }, see Figure [1] We define the Hamiltonians 
Hf{S) and H?(S) as follows. Hf includes all terms d u S u , u £ Xi and all terms c uv S U S V 
such that (u, v) is a vertical edge that has exactly one end-point in Xi. Similarly Hj includes 
all terms d u S u for u S Yj and and all terms c uv S U S V such that (u,v) is a horizontal edge 
that has exactly one end-point in Yj. It is easy to see that Y^iZo H^(S) + Hf(S) = 2H(S"ff 
This implies that there exists an i and b = x or y such that Hf(S op t) > H(S op t)/t or 
H(S pt) — H^(S pt) < (1 — e)H(S opt), where S^pt is a spin configuration with the minimum 
energy. 



i+t 



j+2t 



i+t 



i+2t - 

Fig. 1 . Sets of vertices on horizontal lines X ; and sets of vertices at vertical lines Yj . Drawn are 
also some vertical edges that are part of Hf and horizontal edges that are part of H v - . 



For any i and b — x,y consider a Hamiltonian H^ ubi (S) = H(S) — H%(S). Note that 
this Hamiltonian describes the Ising model defined on a set of t\fn disconnected strips of size 
- x y/n and e^/n disconnected lines of size \Jn. The latter corresponds to sets of edges that 
have both endpoints in Xi or Y\ and thus contains no linear terms. It means that H b sub ^S) 
is invariant under flipping all the spins in any connected component of Xj (if b — x) or Yi (if 
b = y). Besides, as shown in the previous paragraph, there exists a choice of i and b such that 
H b sub JS opt )<(l-e)H(S opt ). 

*We note that a similar decomposition could be obtained for qubits on the square lattice with additional 
diagonal interactions. In such case we would define four Hamiltonians corresponding to lines of vertical, 
horizontal, diagonal-\, and diagonal-/ vertices with the edges incident on these vertices. 
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Let S' opt be a spin configuration that achieves the minimum of H b sub ^S) for some fixed i 
and b. Note that S' opt assigns values to all the vertices of the lattice. Using the symmetry 
of H b ubi (S) mentioned above, one can choose S' t such that H b (S' opt ) < 0. Indeed, if b = x 
then Hf(S' opt ) changes a sign under flipping all the spins in X, while H^ ubi (S) is invariant 
under this flip (the same argument applies to b — y). Therefore, for any i and b one has 
H(S opt ) < H(S' opt ) < H b ub JS' opt ) and for some i and b one has H b subl {S' opt ) < H b ub JS opt ) < 
(1 — e)H(Sopt)- Thus the minimum energy of H b ub { (S) over all i, b, and S approximates 
H(S opt ) within a factor 1 — e. 

It follows that we can get a PTAS by finding the minimum of H b ubi (S) over all choices 
of i = 0, . . . , t — 1, b = x,y, and all spin configurations S. The running time of this PTAS is 
T = § • [Tgtripi^/n, |) + T st „ p (Vn, 1)] • 0(ey/n), where T strip (r, s) is the running time needed 
to find the optimal solution for the Ising spin glass Hamiltonian on a r x s strip. One can 
easily show that 

Lemma 1 There is a dynamic programming algorithm that computes the optimum solution 
for a r x s lattice-blocks using space 0(2 S ) and time 0(r4 s ). 

Proof: Let Bi denote the lattice B restricted to row i or less (B r = B). Let Sbe{ — 1,+1} 
vector of size s and for i € {1, . . . , r}, let V(i, S) be the optimum value of the Hamiltonian 
restricted to Bi such that the variables on the i-th row have assignment S. The dynamic 
program computes and stores V(i,S) for all i and S starting sequentially from i = 1. This 
suffices as the optimum solution for B is exactly ming V(r, S). 

For i = 1, the quantity V(l, S) can be easily computed for each S. Suppose V(i, S) has 
been computed and stored for all S. For an assignment S of row i + l and an assignment S' of 
row i, let Z(i+ 1, S, S') denote the contribution of all terms to the Hamiltonian corresponding 
to vertices on row i + l, and edges in Bi + \ with at least one end point in row i + l. Since the 
assignment S on row i + l can only affect Z, it follows that 

V(i + 1,5) = min(y(i, S') + Z{i + 1, S, S')). (3) 

Since Z(i + 1, S, S') can be computed in time O(s) and V(i, S') are already stored, computing 
V(i + 1, S) for all S takes time 0{s2 s -2 s ). We can speed up this procedure somewhat to 0{A S ) 
by considering the assignments S' in say the Gray code order (where successive assignments 
differ in exactly 1 variable), and hence only O(l) work needs to be done per assignment 5". 

□ 

Thus using the dynamic program leads to a PTAS with an overall running time T = 
0(n4x). 

2.2 Quantum Speed- Up? 

Let us consider how a quantum computer could improve these running times. A possible 
application of quantum searching is inside the dynamic programming. Since many classical 
approximation algorithms rely on dynamic programming, this could be an important area of 
applications for quantum searching. For each row of s spins the dynamic program performs 
2 s minimizations and the minimum is over a function which takes 2 s values. There exists 
a quantum algorithm for finding the minimum of a function [21j which provides a square- 
root speed-up over a brute-force classical minimization. This algorithm uses the Grover 
search algorithm as a subroutine. However in the dynamic program, part of the input to this 
function is stored which implies that this is a problem of searching in a real database. In 
the database setting one has to consider the additional time/hardware overhead in accessing 
the spatially extended database. Optical or classical wave implementations of this type of 
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searching have been considered, see [55]. Whether this application of Grover's algorithm in 
dynamic programming is of genuine interest will depend how its physical implementation 
competes in practice with the capabilities of classical computers. 

An alternative, less efficient, way to divide up the lattice is to remove edges that connect 
subblocks of L x L spins. This leads to a classical approximation algorithm with running 
time 0(2°( 1 / £ >). In this scenario, a direct quantum search algorithm can help in finding the 
minimum energy of each block, but note that this quadratic improvement would still lead to a 
running time which scales as 2°( 1 / c ' which is worse than the dynamic programming method. 

2.3 General Planar Graphs 

Let us now consider the general case of planar graphs. In the lattice case we used the symmetry 
of the lattice to argue that there exists a small subset of edges such that they have a small 
contribution to the Hamiltonian and removing them decomposes the lattice into small disjoint 
blocks. For general planar graphs we cannot use this argument due to the lack of symmetry 
and hence we argue indirectly. We show that the magnitude of the optimal solution is at 
least a constant fraction of the sum of the absolute values of quadratic terms corresponding 
to the edges. This allows us to find a subset of edges with relatively small weight such that 
removing them decomposes the graph into simpler disjoint graphs for which the problem can 
be solved directly. Since the removed edges have small weight adding them back in does not 
increase the Hamiltonian by too much, irrespective of whether they are satisfied or not. 

2-4 Main theorem and Its Proof 

For a planar graph G = (V, E) let W — J2( u v )eE \ c uv\- The key to our PTAS is the following 
result that shows that the value of the optimal assignment scales linearly with W: 
Theorem 1 If G is planar, then H(S op t) < — W/3. 

We begin with a simple property of the optimal solution (the one that minimizes H(S)) 
that holds for an arbitrary graph. Let us write the Hamiltonian as H(S) = Q(S) + L(S), 
where Q(.) is the contribution of the quadratic (2-local) terms and L(.) is the contribution of 
the linear (1-local) terms. 

Claim 1 There exists an optimal solution S opt such that H(S opt ) < and L(S opt ) < 0. 

Proof: It is clear that J^s^iS) — (this is the traceless condition). This implies that 
there must exist a spin configuration with negative energy and thus H(S op t) is negative. The 
second part can be argued by assuming the contrary (L(S , opt ) > 0) and then noting that the 
solution with opposites signs — S'opt is better than S op t itself. □ 

Thus we note that it suffices to show that minsQ(S') < —W/3 since Claim Q] shows that 
mills H(S) < mins Q(S). Hence we consider planar graphs with only quadratic terms in 
H(S), i.e. we assume that d u = for all u. Recall that Bieche et al. [7] has shown that this 
problem can be solved exactly in polynomial time. Our proof of Theorem [T] builds on the 
ideas of Bieche and so we first describe these ideas. We begin with some notation. 

A graph is planar if it can be drawn in the plane such that no edges cross. This drawing 
defines disjoint regions in the plane that are called faces. A cycle in a graph is a collection of 
edges (tti, U2), (u2, Us), . . . , (u^-i, ug) where Ug = u\. Given two cycles C± and C2, their sum 
G\ © C2 is defined as the symmetric difference of C\ and C2 . The faces of planar graph form 
a cycle basis, that is, every cycle can be expressed as a sum of faces. Given an assignment S, 
an edge (u, v) is called unsatisfied if u and v are not assigned according the sign of c uv (i.e. if 
c U v > but S U S V = 1 or if c uv < but S U S V = — 1). A face F is called frustrated if it contains 
an odd number of edges with positive weight. A key observation is that for any assignment, a 
frustrated face must always contain an odd number (hence at least one) of unsatisfied edges. 
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Conversely, if J is a subset of edges such that each frustrated (resp. non-frustrated) face 
contains exactly an odd (resp. even) number of edges in J, then there is an assignment S 
such that the unsatisfied edges are exactly those in J. Thus, Q(S) = —W + 2J2( U v)ej \°uv\- 
Bieche et al. |7, showed that for planar graphs finding such a set J with minimum weight 
is equivalent to finding the minimum weight T-join in the dual graph G* (these terms are 
defined below). Let us remark that the total number of frustrated faces is even because each 
edge (u,v) with c UtV > has exactly two adjacent faces. 

For a planar graph G, its dual graph G* is defined as follows: G* has a vertex for each face 
in G. Vertices u and v in G* are connected by an edge if and only if the faces corresponding 
to u and v in G share a common edge. Given a G, the dual G* is not necessarily unique (it 
depends on the drawing of G), however G** = G. A subset of edges E' is called a cut-set 
if removing them disconnects the graph into two or more components. For planar graphs, 
every cut-set in the dual graph G* corresponds to a cycle in G. We shall see below that sets 
of unsatisfied edges J in G can be identified with T-joins in the dual graph G* , where T is a 
set of frustrated faces. 

Let G = (V, E) be a graph with edge weights, and let T be a subset of vertices T C V 
such that |T| is even. A T-join is a collection of edges J such that each vertex in T is adjacent 
to an odd number of edges in J and each vertex in V \ T is adjacent to an even number of 
edges in J. The minimum weighted T-join problem is to find a T-join with minimum weight, 
and this can be found in polynomial time using matchings. As we mentioned above finding 
the optimal assignment is equivalent to finding the minimum weight T-join in G* where T 
to be the set of vertices corresponding to the frustrated faces in G, and where an edge e 
corresponding to (u,v) in G has weight w(e) = \c uv \. 

We will use a polyhedral description of T-joins. For a subset of edges J, let v( J) denote 
the corresponding |Tv| -dimensional incidence vector (with 1 in the i-th coordinate if edge i 
lies in J and otherwise). For a subset of vertices X, let S(X) denote the set of edges with 
one end point in X and other in V \ X . Given a graph G and the set T, we say that a subset 
of edges J is an upper T-join if some subset J' of J is a T-join for G. Let P be the convex 
hull of all vectors v( J) corresponding to the incidence vector of upper T-joins. P is called the 
up-polyhedra of T-joins. Edmonds and Johnson [23] gave the following exact description of 
P (see the book by Schrijver [21], Chapter 29, par. 1-6 for further details). 

^2 x ( e ) - 1 ' for a11 sets x s - t - \ x n r i is odd ' ( 4 ) 

eeS(W) 

< x(e) < 1 for all edges e. (5) 

This implies that any feasible solution x to the system of inequalities above can be written 
as a convex combination of upper T-joins, i.e. x = ^ aiv ( Ji) where < ck^ < 1 and J^. oti = 1. 
In particular this implies that 

Corollary 1 // all the edge weights w(e) are non-negative, then given any feasible assignment 
x(e) satisfying the inequalities above, there exists a T-join with cost at most ^2 e w(e)x(e). 

We are now ready to prove Theorem [T] 
Lemma 2 Let G be a simple ( with no multiple edges between the same pair of vertices ) planar 
graph with weights \c uv \ on the edges, and let G* be its dual graph. For any subset of vertices T 
of G* such that \T\ is even, the minimum weighted T-join has weight at most Q2( u v) l c ™|)/3- 
Proof: Each cut-set J of G* corresponds to a cycle in G. Since G is simple, each cycle has 
length at least 3, and hence each cut-set J of G* contains at least 3 edges. Consider the 
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assignment x(e) = 1/3. It clearly satisfies Eq. ©. Moreover it also satisfies Eq. (0| as 
5(X) > 3 for all X C V, X ^ V. The result then follows from Corollary Q] 

□ 

Lemma [2] implies Theorem [1] immediately, since the value of the optimal assignment is — W 
plus twice the weight of the optimal T-join which is at most — W + 2W/3 = ~W/3. Observe 
that Theorem Q] is tight, as seen from the example where G is a triangle with edge weights 
+1, +1 and — 1. Here W = 3, but the optimal spin assignment has value — 1. The condition 
that G is simple is necessary. Otherwise, consider the graph on two vertices with two edges, 
one with weight —1 and other with weight +1. Here W — 2, but the optimal assignment has 
value 0. 

2.5 The Approximation Algorithm 

Given Theorem [TJ the PTAS follows using some ideas in [20] • We begin by describing the 
notions of p-outerplanar graphs and tree- widths. An outerface of a planar graph drawn in the 
plane is the set of edges that constitute the boundary of the drawing. For a tree, the outerface 
is the set of all the edges. An outerplanar or 1-outerplanar graph is a planar graph that has 
an embedding in the plane with all vertices appearing on the outerface. Hence examples of 
1-outerplanar graphs are trees or the graph consisting of two cycles that share a common 
vertex. One can define a p-outerplanar graph recursively as (see the example in Figure [2]) : 
Definition 2 (p-outerplanar graphs) A p- outerplanar graph is a planar graph that has an 
embedding in the plane such that removing all the vertices on the outer face gives a (p — 1)- 
outerplanar graph. 



Fig. 2. A planar graph which is drawn as a 4-outerplanar graph. 

The notion of tree decompositions (TD) was introduced by Robertson and Seymour [25] , 
see also [20, 26 . Roughly speaking, a TD allows one to map the Ising spin Hamiltonian H(S) 
on a graph G = (V, E) to a new spin Hamiltonian H tree (Q) that depends on spins ©t € {0, l} b 
living at vertices t of some tree T. The spin O t represents the "opinion" that a vertex t has 
about value of spins S u in some subset of vertices B t C V called a bag. Accordingly, © t 
may take 2 b values, where b is the number of spins in B t . Suppose we can form m bags 
B\ , . . . , B m C V associated with the m vertices of T such that (i) every vertex u G V is 
contained in some bag B t ; (ii) for every edge (u, v) S E, some bag B t contains both u and v; 
(iii) a set of bags containing any given vertex of G forms a subtree of T. These data specify 
a TD of G with a size m and a tree-width 6 — 1. The rules (i),(ii) guarantee that one can 
distribute the terms c uv S u S v and d u S u of the Hamiltonian H(S) over the bags B\, . . . , B m 
such that every term appears in exactly one bag. This distribution defines Hamiltonians 
H t ree,t(@t) = J2(u.v)eB t c U vS u S v + J2ueB t d uS u , where the spins S u are determined by 9*. 
Define H tree (Q) = Y^tLi-B-tree,t(®t)- For every edge (s,t) £ E(T) let us say that O s and 
Ot are consistent on (s, t) iff their opinions about any spin u £ B s D B t agree. The rule 
(iii) guarantees that O is consistent with some spin configuration S iff (O s , Ot) are consistent 
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on every edge of T. Accordingly, ming H (S) — minQ i7 tree (0), where min' means that the 
consistency condition on every edge of T is imposed. The optimal solution can be found very 
efficiently using the standard dynamic programming approach since the problem is defined 
on a tree. It requires a running time of 0(m4 b ), see |26) . 

It is known that a p-outerplanar graph has a TD with a size m = 2n — 1 and a tree- width 
at most 3p — 1. Such a TD can be computed in time 0(pn), see |27j . Summarizing, the 
minimum energy of the Ising spin glass Hamiltonian on a p-outerplanar graph with n vertices 
can be found in time T — 0(n2 6p ). 

Our algorithm works as follows. Given the planar graph G, one first constructs a drawing 
of G in the plane. This can be done in linear time, using for example the algorithm of Hopcroft 
and Tarjan [25]. This gives an outerplanar decomposition of G, see e.g. Fig [2] Say it is h- 
outerplanar (h could be as large as 0(n)). Partition the vertices into levels V\, . . . , Vh where 
Vi is the outer face and Vi is the outer face obtained by removing V\, . . . , VJ_i. Let E t be a 
set of edges that go from Vi to Vi+\, For j = 0, . . . , t— 1, let Gj be the union of sets Ei for all 
i = j (modulo t). (Recall that t = 1/e.) As each edge lies in at most one set Gj, there exists 
some index j such that the sum of \c uv \ over all edges in Gj is at most eC. Remove all the 
edges in Gj from the graph G. This decomposes G into a disjoint collection of i-outerplanar 
graphs Fx, F%, ■ ■ ■ , F f h. We find the minimum energy separately on each of these subgraphs. 

Now consider the quality of the solution obtained for the decomposed problem. Let 
H su b,j{S) = H{S) - J2( u ,v)&G c uvS u S v and let the optimal solution for H sub ^ be S' opt . By 
the reasoning above there exists j such that 

HsubAS'opt) < H subj (S opt ) < H{Sopt) + eW. (6) 

Furthermore, H(S' opt ) < H subd (S' opt )+eW. Thus H (S' opt ) < H(S opt ) + 2eW < (l-6e)H(S opt ) 
by Theorem[T] It follows that we can get a PTAS with a relative error 6e by trying all possible 
j = 0, . . . , t— 1 and choosing the optimal solution S' opt that yields the smallest value of H(S' opt ). 

For a fixed j finding S' opt requires time Tj = ~£ali 0(|F a |2 6t ) = 0(en2 6t ). Thus the 
overall running time of the PTAS is T = £)Lo = 0(n2 6 / e ). Choosing 5 = 6e this implies 
that the algorithm obtains an assignment with the energy at most (1 — S)H(S op t) in time 
0{n2 36 / s ). 

3 Quantum Ising Spin Glass 

The following theorem is the quantum equivalent of Theorem [1] and will be instrumental in 
proving our results: 

Theorem 2 The minimum eigenvalue X(H) of H for a planar graph X(H) < — ^2 U ||£ u ||/5 — 
W/(5 ■ 3 5 ), where W = E (u ,„ )eB \\Q«A- 

Proof: The strategy will be to upper bound X(H) by X sep (H), where A sep is the minimal 
energy achieved on tensor products of states |0), |+), |— ), | + i), \ — i), where 

|±) = (|0)±|1»/V2, |±zH(|0)±z|l»/V2. 

It is enough to prove that X sep is an extensive quantity and this can be achieved using the 
classical result, Theorem[T] Let us first prove the Theorem for the special case then L = 0, that 
is H — Q involves only interactions quadratic in Pauli operators. For every edge (u, v) G E 
the interaction Q u ^ v generally involves all 9 combinations of Pauli operators. We will choose 
one of them that has the largest magnitude and call it a dominating coupling (ties are broken 
arbitrarily). For example, if Q u . v — 3X U ®Y V — 4Z U ®X V , then the edge [u, v) has dominating 
coupling — AZ U ® X v . We have 
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Lemma 3 Suppose Q u , v has a dominating coupling c uv P" (g> Pfi , where P a , P 13 € {X, Y, Z}. 
Then 

\C U v\ > g\\Qu,v\\- ( 7 ) 

Proof: Indeed, otherwise the triangle inequality would imply ||<3u,i;|| < 9|c„„| < ||Q ut; ||. 

□ 

We shall now partition the edges E into several subsets E = UjEj, such that the dominating 
couplings in each subset Ej commute with each other, that is, the sum of dominating couplings 
over Ej is equivalent to a classical Ising Hamiltonian up to a local change of basis. First of all, 
since G is a planar graph we can color its vertices by {1,2,3,4} such that adjacent vertices 
have different colors. A map / : V — > {X, Y, Z} that assigns a Pauli operator to every vertex 
of G will be called a Pauli frame if f(u) depends only on a color of u. Consider the following 
Pauli frames. 



1 


2 


3 


4 


X 


X 


X 


X 


X 


Y 


Z 


Y 


X 


Z 


Y 


Z 


Y 


X 


Z 


Z 


Y 


Y 


Y 


X 


Y 


Z 


X 


Y 


Z 


X 


Y 


Y 


Z 


Y 


X 


Z 


Z 


Z 


z 


X 



This table forms an orthogonal array of strength two with alphabet {X, Y, Z}, that is every 
pair of columns contains every possible combination of two Pauli operators exactly one time. 
Let /i, . . . , fg be the corresponding Pauli frames. Denote Ej a subset of edges (u, v) € E such 
that (u, v) has a dominating coupling 

Then we conclude that 

Ej n E k = if j^k, and u| =1 Ej = E. (8) 

For every Pauli frame fj we can introduce a classical Ising Hamiltonian Qj obtained from 
Q by restricting the whole Hilbert space to classical states in the Pauli frame fj (that is, 
if fj(u) = X for some qubit u, we allow this qubit to be in either of states |+) and |— ); if 
fj(u) = Z, we allow u to be in cither of states |0), |1), e.t.c). By definition, 

A(Q) < X(Qj). (9) 

Note that for every edge (it, v) € Ej the dominating coupling in Q u>v is diagonal in the Pauli 
frame fj. Thus, applying Theorem [T] to Qj we obtain 

KQj)<-\ E l c -l<-^ E IIGim.II> (io) 

where the second inequality follows from Lemma [3] Combining Eqs. (|8I9I10[) we arrive to 

1 9 1 

a(q)<-E a (^)<-^ E HGtM-ii- ( n ) 
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It remains to generalize this bound to the case L ^ 0. We can show that, similar as in the 
classical case (see Claim [I), the following holds: 

Lemma 4 One can choose a ground-state |^>o) of Q such that (^ol-^l^o) < 0. 

Proof: Consider an anti-unitary operator W (known as Kramers- Wannier duality) such that 

w\4>) = ®Y u \p), 

uev 

where the complex conjugation is performed in |0), |1) basis. One can check that P° W — 
— WP" for any single-qubit Pauli operator P". It follows that 

QW = WQ, LW = -WL. 

Thus, if one defines a state |0o) = W \ipo), one gets 

(<h\L\iM = -M)|£|lfo>, (<h\Q\<h) = (ifo\Q\1>o). 

It follows that \cf>o) is also a ground state of Q and one of the expectations values {(po\L\<po) 
or (ipo\L\tpo) is non-positive. 

□ 

Now define 5 states p\ , . . . , p§ such that 

(i) For j = 1, ... ,4 a state pj sets qubits u E V of color j to the ground state of L u ; all other 
qubits are set to the maximally mixed state. 

(ii) p5 is a ground state of Q such that Tr(p 5 L) < 0. 
Then one has the following inequalities 

Tv{Q Pl )=0 for j = l,...,4, and Tr(Qp 5 ) = X(Q). (12) 
Tr(L Pj ) = - II L "H for J = 1, - • - = 4, and Tr(ip 5 )<0. (13) 

u : co\or(u)—j 

Let p be the uniform probabilistic mixture of pi, . . . ,^5. Then Eqs. (|12ll3p imply that 

Tr(pff)<-i]T||L M ||-^ ^ HQ^II. 

uev {u,v)ge 

□ 

3. i T/ie Approximation Algorithm for the Quantum Problem on a Planar Graph 
with Bounded Degree 

We use the following result of Klein, Plotkin and Rao |29j on decomposing planar graphs. 
Theorem 3 Let G = {V, E) be an undirected planar graph with non-negative edge weights, 
and let W denote the total edge weight. Then, given any e > 0, there is a subset of edges 
E' with total weight at most eW such the removing these edges decomposes the graph G into 
components each of which has weak diameter at most c/e with respect to G. Here c is fixed 
constant independent of e. 

Here is the definition of a weak diameter: 
Definition 1 Weak Diameter: Let G = (V, E) be a graph and let G 1 — (V, E 1 ) be a 
subgraph of G. We say that G' has weak diameter d with respect to G, if for any two vertices 
V, W G V , their distance in G is at most d. 
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If G is a planar graph with maximum degree d, Theorem[3]implies that each component can 
have at most dP^ 1 /^ vertices and hence the minimum eigenvalue problem for a Hamiltononian 
H restricted to every component can be solved in time 2 d ° l / ' . 

There is a linear time algorithm to determine the set of edges E' . The algorithm works as 
follows: Let 5 = e/3. Root the graph G at arbitrary vertex and consider the breadth first tree. 
A vertex is said to be at level i, if it is at distance i from the root. For j = 0, . . . , 1/5 — 1, let 
Ej denote the set of edges that connect two vertices at level i and i + 1 where i = j (modulo 
1/(5). Choose the set Ej with least weight and remove these edges from G. Let G\ denote 
the obtained graph. Now consider each of the components of G± and apply this procedure 
again to obtain the graph G2. Finally, apply the same procedure to G2 to obtain G3 (one 
applies the procedure three times because planar graphs K3 3 minor free). The result of Klein, 
Plotkin and Rao [S3] shows that G3 has weak diameter at most 0(l/<5) = 0(l/e). Moreover 
the weight of edges removed is at most 35 = e fraction of the total weight. 

Let H = J2(u v)eE Qu.v + Yluev ^ e a <l uan t um Ising spin glass Hamiltonian defined on 
a graph G — (V 7 E). Define a weight associated with an edge (u, v) as ||Q„.t,||. Let H sub be 
a Hamiltonian obtained from H by keeping all the linear terms L u and the quadratic terms 
Q u ,v associated with edges of a subgraph G3 defined above. By definition of G3 one has 
\\H - H sub \\ < eW and thus \X(H sub ) ~ X(H)\ < eW. Theorem H implies that \X(H sub ) - 
X{H)\ < ce|A(iJ)| for some numeric constant c. Thus one can approximate X(H) with any 
fixed relative error e in time n2 d ( 1 ' . 

4 Quantum Ising Spin Glass Problem on a Star Graph 

The construction of PTAS for classical Hamiltonians on planar graphs presented in Section [5] 
relies on the fact that the classical problem on a tree (or any graph of constant treewidth) can 
be solved efficiently using the dynamic programming. Unfortunately, this method does not 
work for quantum Hamiltonians. In this section we develop a new technique that allows one 
to obtain a PTAS for the quantum Ising spin glass Hamiltonian on a star graph — a tree that 
consists of n+l vertices with one vertex having degree n and n vertices having degree 1. The 
corresponding graph is G = (V, E), where V = {0, 1, . . . , n} and E — {(0, 1), (0, 2), . . . , (0, n)}. 
We shall refer to spins sitting at vertices 1, 2, . . . , n as bath spins and the spin sitting at the 
vertex as central spin (by analogy with the central spin problem studied in condensed matter 
physics [32] )■ Let the Hamiltonian be 

n 

H = H + Y J Hoj (14) 
i=i 

where Hq is a linear term acting on the central spin and Hqj is the interaction between the 
central spin and j-th bath spin (which can include both quadratic and linear terms) . 

Theorem 3 Suppose there exist constants < a < b such that a < ||7?o.j|| < b for all 

j. Then one can approximate the smallest eigenvalue X(H) with a relative error e in time 

,-o(i) 

71 

In the rest of this section we prove Theorem [3] We start from proving that X(H) can be 
computed exactly in time poly(n) as long as the number of distinct interactions Ho.j is bounded 
by a constant. We shall use 

Lemma 5 Suppose the interactions Hqj are the same for some subset of bath spins S C 
{1, . . . , n}. Then one can choose a pure ground state of H symmetric under permutations of 
spins in S. 
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Proof: Without loss of generality S — {1,2,..., fe}. Let |^) be a ground state of H. 
Denote Wj k the swap of qubits j and k. Assume that \^>) is not symmetric under permutations 
of spins in S. Without loss of generality, W\z \ *&) ^ |^)- There are two cases: (i) W1.2 \ *f?) 
is proportional to |^). Then W\$ \^) — — since W\$ has eigenvalues ±1. Therefore 
I*) = |*~)i2 ® |* e ; se ), where |<f~) = ^(|01) - 1 10) ) is the singlet state and |* e ; se ) is some 
state of spins {0, 3,4,..., n}. The second case is (ii) W\ t 2 I*) and |\&) are linearly independent. 
Then the anti-symmetrized state \*&') = \ty) — W12 \^) is non-zero. On the other hand, |\&') is 
a ground state of H since W12 commutes with H . We conclude that |W) = | \[/ } 12 ® Y^eise)- 
In both case we conclude that H has a ground state |\&) = | \I/ ) 12 ® Y&eise)- 

Since the energy of a state depends only upon the reduced density matrices poj, we can 
replace the antisymmetric singlet I* - ) by the symmetric EPR state \^ + ) without changing 
the energy. On the other hand, any state with energy equal to the ground-state energy must 
be a ground state. We conclude that H has a ground-state |^) = |4' + )i2 <S> \^else)- 

By iterating the arguments above one concludes that H has a ground state 

I*) = l*ia) ® • ® l*fc-l,fc> ® = |*S) ® l^else)- (15) 

where l^eZse) is some state of all spins j <£ S (if k is odd then there will be one unpaired spins 
in S; this will not change the arguments below very much). Now we can symmetrize |^) by 
brute force method. Let 

rGSfc 

be the projector onto the symmetric subspace, where W(t) is the unitary operator implement- 
ing a permutation r of k spins in S. Note that the state \^s) in Eq. (fT5|) has non-negative 
amplitudes in the standard basis. Therefore W(t) \^s) also has non-negative amplitudes. 
Therefore II \&s) ^ 0, and, accordingly, \^>') — II5 eg) I e i se |W) is a non-zero state symmetric 
under permutations of spins in S. On the other hand, since W(t) commutes with H, \ty') is 
a ground state of H. 

□ 

This result implies that we can look for a ground state that "occupies" only a \S\ + 1 
dimensional subspace of the 2' s ' dimensional Hilbert space describing spins in S. If we have 
M distinct interactions, the dimension of the space in which the optimization takes place is 
n^f 1 (|S'i| + 1) < n M which is polynomial in n. Thus the optimization problem for constant 
M can be solved efficiently in n. 

In order to map the general problem onto one in which we have constant number of distinct 
interaction, we apply a coarse-graining procedure to the general Hamiltonian Eq. (|14j) . One 
can show 

Lemma 6 For any < a < 1 define a set M a of 2-qubit Hamiltonians H satisfying a < 
\\H\\ < 1. For any e > there exist 2-qubit Hamiltonians G%, . .., Gm, M = 0((ae)~ 32 ) such 
that min a \\G a - H \\ < e \\H\\ for any H G M a . 

Proof: It is enough to satisfy min Q \\G a — H\\ < ea. A 2-qubit Hamiltonian satisfying 
\\H\\ < 1 lives ina2x2x---x2 cube in R 32 . Construct ea-mesh, count the number of points. 

□ 

Now we are ready to prove Theorem |3] Without loss of generality, b — 1 (otherwise 
multiply H by 6 _1 ). Applying Lemma [5] to every interaction Hqj one can partition the 
n bath spins into M = 0((ae)~ 32 ) subsets S\,...,Sm such that \\H j — G a \\ < e\\H _j\\ 
for all j E S a - We define a coarse-grained Hamiltonian H = Hq + ^ a=1 J2jes G a \j]- 
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Here the notation G a [j] means that G a acts on the spins and j. We have \\H — H\\ < 
e^" =1 ||-ffo,j|| < e5-3 5 |A(ff)| where the second inequality follows from Theorem^ Therefore, 
\X(H) — X(H)\ < ce\X(H)\ for some numeric constant c. The classical PTAS will find the 
ground-state (that is, a poly(n)-sized classical description of this state) and the ground-state 
energy X(H) of the coarse-grained Hamiltonian. Lemma[5]implies that it requires time rf ° ( . 

5 Discussion 

An important open question is whether there exists a classical or quantum PTAS for the gen- 
eral quantum Ising spin glass problem on planar graphs. It is clear that some new techniques 
will be needed to settle this problem. A simpler problem in this realm would the quantum 
Ising spin glass problem on a tree with unbounded degree. Note that even in the simplest 
case of star graphs the existence of PTAS for the quantum problem is not proved (Theorem [3] 
assumes the constant lower and upper bounds on the norm of interactions ifo.j)- 

One interesting approach to address these quantum problem may be to consider quantum 
or classical algorithms that output the thermal state e~ H l T jZ at temperature T for the Ising 
spin glass problem. Such thermal state will typically provide a PTAS for the ground-state 
energy problem. One can prove this by showing that the average energy (H )t = Tri? e~ H ^ T Z 
is bounded as \ (H)t — X(H)\ < 2nT. This bound follows from the fact that for the free-energy 
F(T) we have \F(T) - (H) T \ < nT and \F(T) - X(H)\ = \F(T) - F(0)| < nT. When the 
ground-state energy X{H) scales with n (e.g. for bounded weights \c uv \ < c, HQ^H < c one 
gets this from Theorems Q] and [5]) , the error in the approximation can be made eA(-ff) for T = 
0(e). This also shows that finite but small temperature implementation of adiabatic quantum 
computation will generally provide a PTAS-approximation to the ground-state energy problem 
(assuming that, say, for bounded-degree graphs beyond the planar ones, the ground-state 
energy will be extensive, scales with n). 

The difference between classical and quantum behavior on tree graphs is also witnessed by 
the fact that the algorithm of classical belief propagation (for zero temperature this essentially 
corresponds to dynamic programming) converges efficiently on trees, whereas quantum belief 
propagation will only work when additional conditions are fulfilled [301 131] . It is expected 
that for bounded-degree trees the quantum belief propagation algorithm of [31] at finite 
temperature T will give rise to a PTAS. 

Acknowledgements 

SB and BMT acknowledge support by NSA and ARDA through ARO contract number 
W911NF-04-C-0098. 

1. K.H. Fischer and J. A. Hertz. Spin Glasses. Cambridge University Press, Cambridge, 1991. 

2. Y. Fu and P. W. Anderson. Application of statistical mechanics to NP-complete problems in 
combinatorial optimisation. J. Phys. A: Math, and Gen., 19(9), pp. 1605-1620 (1986). 

3. S. Kirkpatrick and B. Selman. Statistical physics and computational complexity. 
http : //www. cs . Cornell . edu/selman/papers/pdf /01 . condensed. physics .pdf 

4. F. Barahona. On the computational complexity of Ising spin glass models. J. Phys. A: Math, and 
Gen., 15, pp. 3241-3253 (1982). 

5. R.I. Oliveira and B.M. Terhal. The complexity of quantum spin systems on a two-dimensional 
square lattice, http://arxiv.org/abs/quant-ph/0504050 

6. D. Aharonov, D. Gottesman, S. Irani, and J. Kempe. The power of quantum systems on a line. 
In Proceedings of 48th FOCS (2007) , [http : //arxiv . org/abs/0705 . 4077 1 

7. L. Bieche, R. Maynard, R. Rammal, and J. P. Uhry. On the ground states of the frustration model of 
a spin glass by a matching method of graph theory. J. Phys. A: Math, and Gen., 13, pp. 2553-2576 



15 



(1980). 

8. E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser. Quantum computation by adiabatic evolution. 
http : //arxiv . org/abs/quant-ph/ 0001 106 

9. A. Ambainis and O. Regev. An elementary proof of the quantum adiabatic theorem. 
http : //arxiv. org/abs/quant-ph/041 1152 

10. W. van Dam, M. Mosca, and U. Vazirani. How powerful is adiabatic quantum computation? In 
Proc. of FOCS, pp. 279-287 (2001). 

11. E. Farhi, J. Goldstone, S. Gutmann, and D. Nagaj. How to make the quantum adiabatic algorithm 
fail http : / / arxiv . org/ abs/quant-ph/ 0512159 

12. J. Brooke, D. Bitko, T. F. Rosenbaum, and G. Aeppli Quantum Annealing of a Disordered Magnet 
Science 284 (5415), 779 (1999) DOI: 10.1126/science. 284.5415. 779 

13. R.D. Somma and S. Boixo and H. Barnum and E. Knill Quantum Simulations of Classical An- 
nealing Processes http : //www. citebase . org/abstract?id=oai : arXiv. org: 0804 .1571 

14. G.Rose, http://dwave.wordpress.com 

15. V.V. Vazirani. Approximation Algorithms. Springer- Verlag, Berlin, Germany (2001). 

16. P. Alimonti and V. Kann. Some APX- completeness results for cubic graphs. Theor. Comput. Sci., 
237, pp. 1-2 (2000). 

17. M. Jerrum and A. Sinclair. Polynomial-time approximation algorithms for ising model (extended 
abstract). In Automata, Languages and Programming, pp. 462-475 (1990). 

18. R. J. Lipton and R. E. Tarjan. Applications of a planar separator theorem. SIAM Jour, on Comp., 
9, pp. 615-627 (1980). 

19. B. S. Baker. Approximation algorithms for NP-complete problems on planar graphs. Journal of 
the ACM (JACM), 41, pp. 153-180 (1994). 

20. S. Khanna and R. Motwani. Towards a syntactic characterization of PTAS. In In Proc. of STOC, 
pp. 329-337 (1996). 

21. C. Diirr and P. H0yer. A quantum algorithm for finding the minimum. 
|http : //arxiv ■"or g/abs/quant-ph/9607014] 

22. P.G. Kwiat, J.P. Mitchell, P.D.D. Schwindt, and A.G. White. Grover's search algorithm: An 
optical approach. J. Mod. Opt., 47, pp. 257-266 (2000). 

23. J. Edmonds and E. Johnson. Matchings, Euler tours and the Chinese postman problem. Math. 
Programming, 5, pp. 88-124 (1973). 

24. A. Schrijver. Combinatorial Optimization: Polyhedra and Efficiency. Springer- Verlag (2003). 

25. N. Robertson and P. D. Seymour. Graph minors II: Algorithmic aspects of treewidth. Journal of 
Algorithms, 7, pp. 309-322 (1986). 

26. H. L. Bodlaender. A tourist guide through treewidth. Acta Cybernetica, 11, pp. 1-22 (1993). 

27. H.L. Bodlaender. A partial k-arboretum of graphs with bounded treewidth. Theor. Comp. Sci., 
209, pp. 1-45 (1998). 

28. J. Hopcroft and R.E. Tarjan. Efficient planarity testing. Jour, of the ACM, 21, pp. 549-568 
(1974). 

29. Philip Klein, Serge A. Plotkin, and Satish Rao. Excluded minors, network decomposition, and 
multicommodity flow. In Proceedings of STOC, pp. 682-690 (1993). 

30. M. Leifer and D. Poulin. Quantum Graphical Models and Belief Propagation. 
http T//arxiv.org /abs/070 8. 1337| 

31. M. B. Hastings. Quantum belief propagation. Phys. Rev. B RAPIDS, 76, p. 201102 (2007). 

32. M. Gaudin. Diagonalisation d'une classe d'hamiltoniens de spin J. de Physique 37, p. 1087 
(1976). 



10 



